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ABSTRACT 

Ultra-strong emission-line galaxies (USELs) with extremely high equivalent widths (EW(H/3) > 
30 A) can be used to pick out galaxies of extremely low metallicity in the z = — 1 redshift range. 
Large numbers of these objects are easily detected in deep narrow band searches and, since most 
have detectable [OIIIJA4363, their metallicities can be determined using the direct method. These 
large samples hold out the possibility for determining whether there is a metallicity floor for the 
galaxy population. In this, the second of our paper s on t he topic, we describe the results of an 
extensive spectroscopic follow-up of the iKakazu et al.l (|2007l ) catalog of 542 USELs carried out with 
the DEIMOS spectrograph on Keck. We have obtained high S/N spectra of 348 galaxies. The two 
lowest metallicity galaxies in our sample have 12-|-log(0/H) = 6.97±0.17 and 7.25±0.03 - values 
comparable to the lowest metallicity galaxies found to date. We determine an empirical relation 
between metallicity and the R23 parameter for our sample, and we compare to this to the relationship 
for low redshift galaxies. The determined metallicity-luminosity relation for this sample is compared 
with that of magnitude selected samples in the same redshift range. The emission line selected galaxies 
show a metal-luminosity relation where the metallicity decreases with luminosity and they appear to 
define the lower bound of the galaxy metallicity distribution at a given continuum luminosity. We also 
compute the Ha luminosity function of the USELs as a function of redshift and use this to compute 
an upper bound on the Lya emitter luminosity function over the z = — 1 redshift range. 
Subject headings: cosmology: observations — galaxies: distances and redshifts — galaxies: abundances 
— galaxies: evolution — galaxies: starburst 



1. INTRODUCTION 

Developing a large sample of low metallicity galax- 
ies is of considerable interest for clues it can provide 
to the early stages of galaxy formation and chemical 
enrichment — such as whether forming galaxies have 
a baseline metallicity that reflects the early chemical 
enrichment of the intergalactic medium. The most 
metal-poor systems currently known are the low red- 
shift blue compact emission-line galaxies such as I Zw 
18 and SBS 335-052W with measured 12-Hok (0/H) 
of ~ 7.1 - 7.2 (iSargent fc Searia fTOTnl: iThuan fc IzotoTl 
120051: llzotov et al.ll2005D . However, despite enormous ef- 
forts, only a few dozen xMPGs (objects with (12-|-log 
(O/H) < 7.65 or Z < Zq/12; [Kniazev et al. 2003: 
Izotov et al.l l2Q06aD are known fe.g.. IOevll2006t llzotovl 



2006bH . There are too few of these for detailed studies of 



the metallicity distribution function. 

In the first paper of the present series we showed that 
large samples of low metallicity objects can be found by 
searching for ultra-strong emission-line galaxies (USELs) 
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using very deep narrowband surveys (iKakazu et al.l 
l2007h . The narrowband method has many advantages; 
it is an efficient way to develop very large samples of ob- 
jects and, since it probes to much deeper line-flux limits 
than objective prism or continuum surveys, we can study 
populations out to near redshift z 1 where the cosmic 
star formation rates peak. 

In the present work we report on the detailed spec- 
troscopi c follow-up of the ca talog of 542 USEL g alaxies 
given in IKakazu et al.l (|2007f ). The IKakazu et aiT sample 
was chosen from a set of narrowband images obtained 
with the SuprimeCam mosaic CCD camera on the Sub- 
aru 8.2-m using two ~120 A (FWHM) filters centered at 
nominal wavelengths of 8150 A and 9140 A in regions of 
low sky background between the OH bands. The total 
covered area is 0.5 square degrees. The catalog consists of 
all objects with narrow band magnitudes less than 25 and 
a narrow band excess of 0.8 magnitudes for the 8150 A 
filter and of 1.0 magnitudes in the 9140 A filter compared 
to their respective reference continuum bands. The sam- 
ple should contain all galaxies in the fields with observed 
frame equivalent widths significantly greater than 120A 
and 180A in the two filters, and with line fluxes above 
^ 1.5 X 10^^'' erg cm"^ s~^. 

Spectroscopic redshifts have now been obtained for 
299 of these galaxies usi ng multi-obj ect masks with the 
DEIMOS spectrograph (jFaber et al.ll2003D on the 10-m 
Keck II telescope. The observations, flux calibration, 
and equivalent width measurements are discussed in §2, 
where the resulting measured line ratios are also de- 
scribed. In §3 we analyse the metallicities and describe 
an empirical relation for the metallicity as a function of 
the R23 parameter, which we compare with low redshift 
determinations of the relation. The distribution of metal- 
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Fig. 1. — Redshift versus narrow band magnitude in the selection 
filter. The lower histogram (blue line) shows the fraction of objects 
in the magnitude range which have either not been observed or not 
been identified. 

licities as a function of the properties of the galaxies is 
given in §4 where we also compare the results with those 
in magnitude selected samples at the same redshift. In §5 
we show that the results can be used to place upper limits 
on the Lyman alpha emitter (LAE) luminosity function 
in the z = — 1 redshift range, and c ompare this both 
with the local GALEX measurement of 'Dehar veng et aTl 
(|2008 ) and with higher redshift LAE functions. A final 
summary discussion is given in §6. We use a standard 
Ho = 70 km s-i Mpc~\ f},„ = 0.3, Oa = 0.7 cosmology 
throughout the paper. 

2. SPECTROSOPIC OBSERVATIONS 

(|Kakazu et al.l l2007| ) reported spectroscopic observa- 
tions for 161 USELs from their sample. These spec- 
tra were obtained using the Deep Extraga lactic Imag- 
ing M ulti-Object Spectrograph (DEIMOS; iFaber et all 
|2003| ) on the Keck II 10-m telescope in a series of runs 
between 2003 and 2006. Over the 2007 and 2008 period 
we have roughly doubled this spectroscopic sample and 
have also significantly deepened the spectra of many of 
the objects that had previously been observed. 

In order to provide the widest possible wavelength cov- 
erage the observations were primarily made with the 
ZD600 £/mm grating. We used 1" wide slitlets which 
in this configuration give a resolution is 4.5 A, sufficient 
to distinguish the [OIIJA3727 doublet structure. This 
allows us to easily identify [O IIJA3727 emitters where of- 
ten the [Oil] A3 72 7 doublet is the only emission feature. 
The spectra cover a wavelength range of approximately 
5000 A with an average central wavelength of 7200 A, 
though the exact wavelength range for each spectrum 
depends on the slit position with respect to the center 
of the mask along the dispersion direction. The observa- 
tions were not generally taken at the parallactic angle, 
since this was determined by the mask orientation, so 
considerable care must be taken in measuring line fluxes, 
as we discuss below. Each ~ 1 hr exposure was broken 
into three subsets, with the objects stepped along the slit 
by 1.5" in each direction. Some USELs were observed 
multiple times, resulting in total exposure times of up 
to 10 hours. The two-dimensional s pectra were r e duced 
following the procedure described in lCowie et al.l ()1996f ) 
and the final one-dimensional spectra were extracted us- 



FlG. 2. — Redshift distribution of the spectroscopically observed 
sample. The peaks correspond to the positions at which the strong 
emission lines cross the two filters. Thus there are galaxies at 
z = 0.6 and z = 0.8 corresponding to the [OIII]A5007 line lying in 
the 8150 A and 9140 A filters respectively. 

ing a profile weighting based on the strongest emission 
line in the spectrum. 

Our primary goal was to obtain a nearly complete 
set of spectra for the 200 objects in the catalog with 
narrow-band magnitudes brighter than 24, since, for 
these brighter objects, we can obtain high quality spectra 
suitable for detailed analysis. However, we also observed 
a large number of the galaxies with magnitudes in the 
24-25 range. We have now observed 171 of the galax- 
ies brighter than 24. Nearly all of the bright emission 
line candidates which were observed were identified (164 
of the 171 objects). However, three of the objects in 
the sample are stars where the absorption line structure 
mimics emission in the band. The redshift distribution of 
the sample and the fraction of objects which have so far 
been identified is shown as a function of the narrow-band 
magnitude in Figure [T] 

The narrow-band emission-line selection produces a 
mixture of objects corresponding to Ha, [OIII]A5007, 
and [OIIJA3727 and, at the faintest magnitudes (> 24), 
high redshift (z > 5) Lya emitters. There are almost 
no high redshift Lyman alpha emitters at magnitudes 
brighter than 24. The distribution of the redshifts for 
objects other than the LAEs is shown in Figure [H The 
largest fraction of objects corresponds to cases where the 
[OIII]A5007 line lies in the narrow band filters. For the 
present work only the objects selected in [O III] or Ha 
are of interest, since we cannot determine the metallic- 
ities or the Balmer line strengths of the [O II] selected 
samples without futher near infrared spectroscopy. Our 
final sample of galaxies therefore consists of 214 galax- 
ies with redshifts between zero and one of which 189 are 
chosen with [OIII] and 25 with Ha. 125 of these have 
narrow band magnitudes brighter than AB=24. 

3. FLUX CALIBRATIONS 

Generally our spectra were not obtained at the par- 
allactic angle since this is determined by the DEIMOS 
mask orientation required to maximize object placement 
in slits over the entire mask field. Therefore, flux cal- 
ibration using standard stars is problematic because of 
atmospheric refraction effects, and special care must be 
taken for the flux calibra tion. A much more extensive 
discussion can be found in iKakazu et al.l (|2007f) but here 
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Fig. 3. — Portions of the spectrum of the the lowest metaUicity galaxy in the sample showing detected emission lines. This is object 29 
in the IKakazu et aTl I I2007I) catalog of objects selected in the NB912 filter. It lies at a redshift of 0.3931. The emission line features are 
labeled and marked with the solid lines. In panels (a) and (b) we have increased the scale of the vertical axis to show the extremely strong 
H«, [OIII], and H/3 lines. The [Oil] line is only marginally detected in this spectrum and the two [Nil] lines in panel (a) are not seen. 



we focus only on measurements of the relative line fluxes 
from the spectra which we require for the metaUicity 
measurements . 

The spectra were initially relatively calibrated using 
the measured instrument response. Portions of the spec- 
tra of the two lowest metaUicity galaxies in the sample 
are shown in Figures [3] and [H Relative line fluxes can 
be robustly measured from the spectra without flux cal- 
ibration as long as we restrict the line measurements to 
short wavelength ranges where the DEIMOS response 
is essentially constant. For example, one can assume 
that the responses of neighboring lines (e.g. [O IIIJA4959 
and [OIIIJA5007) are the same and therefore one can 
measure the flux ratio without calibration. The present 
problem is considerably simplified by the presence of the 
line near [OIII]A4959 and [OIII]A5007 (Figure [Sljb) 
and Figure m^a)) and the II7 line near [OIIIJA4363 (Fig- 
ure [3]Jc) and Figure HKb)). For these lines, and also for 
the [N II] A6584 line (Figure[3l^a)) we can therefore use the 
Balmer line ratios to provide extinction corrected flux ra- 
tios for the metal lines. To calibrate these lines, we used 
neighboring Balmer lines with the assumption of Case B 
recombination conditions. We cannot so easily do this 
near [O IIJA3727 where the Balmer lines are weak and in 
some cases contaminated (Figure [31[d) and Figure [H^c)). 
Fortunately the [OIIJA3727 is generally extremely weak 
in the spectra and the uncertainty in the calibration has 
little effect on the metal determination. 

For each spectrum we fltted a standard set of lines. 
For the stronger lines we used a full Gaussian fit together 
with a linear fit to the continuum baseline. For weaker 



lines we held the full width constant using the value mea- 
sured in the stronger lines and set the central wavelength 
to the nominal redshifted value. The [O IIJA3727 line was 
fitted with two Gaussians with the appropriate wave- 
length separation. We also measured the noise as a func- 
tion of wavelength by fitting to random positions in the 
spectrum and computing the dispersion in the results. 

The [NII]A6584/Ha, [O III]A5007/H/3, and 
[O II]A3727/[0 IIIJA5007 ratios are shown as a function 
of the narrow band magnitude in the selection filter in 
Figure O The population is relatively uniform in its 
line properties. Nearly all of the galaxies have weak 
[NII]/Ha, (a median ratio of 0.02), weak [OII]A3727 
relative to [OIII]A5007, (a median ratio of 0.22), and 
strong [OIII]A5007/H/3, (a median ratio of 5.2). There 
seems to be little difference between the Ha selected 
population (shown with purple diamonds in Figure [5b) 
and the [O III] selected population (shown with black 
squares). The very weak [Nil] lines, in combination 
with the very strong [O III] , show that these galaxies 
are not excited by active galactic nuclei, and that the 
galaxies have very high ionization parameters and low 
metaUicity, as we shall quantify in the next section. 

4. GALAXY METALLICITIES 

The spectra are of variable quality, reflecting the range 
of magnitudes and exposure times, and, in order to mea- 
sure the metallicities, we need very high signal to noise 
observations. It is also important that the Balmer lines 
are well detected since our flux calibrations rely on the 
neighboring Balmer lines. We therefore restrict ourselves 
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Fig. 4. — Portions of the spectrum of the second lowest metaUic- 
ity galaxy in the sample showing det ected emission lines. This is 
object 269 in the lKakazu et aLl I I2007I) catalog of objects selected in 
the NB912 filter. It lies at a redshift of 0.8175. In (a) we have in- 
creased the scale of the vertical axis to show the extremely strong 
[OIII] and H/3 lines. The emission line features are labeled and 
marked with the solid lines. 



to emitters whose H7 line fluxes are detected above the 
ten sigma level. Among the 25 Ha selected emitters in 
our total spectroscopic sample there are 8 sources with 
spectra of sufficient quality for the metallicity analysis, 
and among the 189 sources selected using [OIIIJA5007 
there are 23 such sources, giving a total sample of 31 
sources for the metallicity analysis. 

As we illustrate in Figure [H nearly all of these sources 
are detected in the [OIIIJA4363 auroral line. Of the 31 
sources, 23 are detected above the 3 sigma level and 9 
above the 5 sigma level. The median value of the ratio 
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Fig. 5.— Flux ratios of [N II] A6584/Ha (panel a), 
[OIII]A5007/H/3 (b), and [O II]A3727/[0 III] A5007 (c) computed 
from the spectra. The line ratios are only plotted for spectra 
where the signal to noise of the H/3 line is greater than five. In 
a small number of cases the line fluxes of the weak [NII]A6584 and 
[OII]A3727 lines in the spectra scatter to negative values. The 
data are plotted against the narrow band magnitude and the error 
bars are one sigma. Ha selected spectra are shown with purple di- 
amonds in panel (b). The red solid curves show the median values 
in the sample. 

of the [OIII]A4363 auroral line to the [OIII]A5007 line is 
0.018 for the 31 objects and there is no indication that the 
values seen in spectra selected with the Ha line, which 
are shown by the purple diamonds in Figure [SI are any 
different from those selected with [O HI] A5007 which are 
shown by the black squares. 

The presence of [OHIJA4363, immediately suggests 
that these are metal-deficient systems but, more im- 
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Fig. 6.— Ratio of the [O III]A4363/[0 III]A5007 fluxes. The 
data are plotted against the narrow band magnitude and the error 
bars are one sigma. Only spectra where the signal to noise of the 
H7 line is above ten are included. The values for the Ha selected 
spectra are shown with purple diamonds and for the [O III] A5007 
selected spectra with black squares. The red solid line show the 
median value in the combined sample of 31 galaxies. The black 
solid line shows the zero ratio. Only three of the sources do not 
have [OIII]A4363 detected above the one sigma level. 



portantly, it allows us to determine the electron tem- 
perature from the ratio of the [OIIIJA4363 line to 
[OIII]AA5007,4959. This procedure is oft en referred to 
as the 'direct' method or Te method (e.g., ISeatonlll975l: 
Pagel et al.lll992t iPilvugin fc ThuanI 120051 : llzotov et all 



To derive Te [O II I] and t he oxygen abundances, 



we used the llzotov et alJ (|2006cl ) formulae, which were 
developed with the lates t atomic da t a and photoioniza- 
tion models. Using the iPagel et al.l ()1992f) calibrations 
with the Te[OII]-Te[OIII] relations derived bv [Garneti3 
(11992), gives consistent abundances within 0.1 dex. The 
[SII]AA6717, 6731 lines that are usually used for the de- 
termination of the electron number density, are beyond 
the Keck/DEIMOS wavelength coverage for our [OIII] 
emitters. Therefore we assumed Ug = 100 cm~'^. How- 
ever the choice of electron density has little effect as elec- 
tron temperature is insensitive to the electron density; 
indeed we get the same results even when we use Ug = 
1,000 cm-3. 

The spectra of the three lowest redshift Ha emission- 
line selected galaxies do not cover the [Oil] A3 72 7 line, 
leaving us with a final sample of 28 galaxies for our 
analysis. Seven of these objects satisfy the definition 
of XMPGs [12 + log(0/H) < 7.65; iKunth fc OstlhH 
[2000]. The lowest metallicity galaxes in the sample are 
KHC912-29 with 12 -I- log(0/H) = 6.97 ± 0.17, and 
KHC912-269 which has 12 + log(0/H) = 7.25 ± 0.03. 
The spectra of these two galaxies are shown in Figures [3] 
and m Their metallicities are comparable to the cur- 
rently known most metal-poor galaxies [I Zw 18 and 
SBS0335-052W; 12 -I- log(0/H) - 7.1 - 7.2]. 

5. DISCUSSION 

5.1. Metallicity and the ionization parameter 

Figure [7] shows the electron temperature sen- 
sitive line ratio, [OIII](A4363)/[OIII]A5007 versus 
[OII]A3727/[OIII]A5007. If we have an estimate of 
the metallicity, as in the present case, we can use the 



O 

o 
in 




0.0 0.1 0,2 0.3 0.4 
[OII]3727/[OIII]5007 



0.5 



Fig. 7.— [OIII]A4959-(-A5007/[O III]A4363 versus 

[OII]A3727/[OIII]A5007 for the 28 galaxies in the metaUic- 
ity sample. One sigma error bars are shown for both ratios. 



[OII]A3727/[OIII]A5007 ratio to estimate the ionization 
parameter q, defined as the number of hydrogen ioniz- 
ing photons passing through a unit area per second per 
unit hydrogen number density. We can see from Figure [7] 
that there is a strong inverse correlation of [OIII]A4363 
and [OII]A3727. Systems with weak [OII]A3727 gener- 
ally have strong [OIII]A4363. 

We have computed the ionization parameters (g) for 
the sample using the parameterized forms of the depen- 
dence oL[OIIlA3727/igHI]A50G7 on q and 12+log(0/H) 
from iKobulnickv fc KewlevI (|2004D . In Figure (Hfa) 
we show the dependence of the q parameter on 
[OII]A3727/[OIII]A5007 and in Figure we show the 
dependence of 12-hlog(0 /H) on [O II] A3727/ [O HI] A4363. 

Their is a clear dependence between the metallicity 
and [OII]A3727/[OIII]A5007. Objects with low metal- 
licity have low [OII]A3727/[OIII]A5007 and all of the 
XlVIPGs have values of this ratio below 0.12. Conversely, 
seven of the thirteen galaxies satisfying this condition 
are XlVIPGs. This is potentially a very valuable tool for 
optimizing the search for the lowest metallicity galax- 
ies at these redshifts since we can focus our spectrosopic 
foUowup on galaxies with weak [OII]A3727 lines. 

The ionization parameter also increases with decreas- 
ing [OII]A3727/[OIII]A5007 (Figure Eta)) though the 
ionization parameters lie in a relatively narrow range for 
log (q)^ 7. 8-8. 5. However, there is a c onsiderable vari- 
ation in the t heoretical models (e.g., llVIcGaugl] Il991l : 
iPilvuginI 120001 ) and the exact shape and normaliza- 
tion are somewhat uncertain. The range of ionization 
parameters is similar to the values found in magni- 
tude limited samples of galaxies at the same redshift 
(jCowie and Bargeij[2008[ ). 

In Figure [9] we show the ionization parameter q as a 
function of the metallicity. It is clear that the lower 
metallicity galaxies have higher ionization parameters. 
However, we again note the uncertainties in the theoret- 
ical modelling. 

5.2. Metallicity and the R23 parameter 

The R23 ratio [/([OIII] A4959) + /([OIII] A5007) -f 
/([OII] A3727)]//(H/3) of Pagel et al. (1979) is one of 
the most frequently used metallicity diagnostics. As is 
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Fig. 8. — (a) The ionization parameter q versus 
[OII]A3727/[Oin]A5007 for the metallicity sample. (b) 12 
+log{0/H) versus [O II]A3727/[0 III] A5007 for the metaUicity sam- 
ple. Only the 25 olajects with [OIII]A4363 detected above the one 
sigma level are shown and one sigma error bars are given in both 
cases. 
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Fig. 9. — The q parameter versus oxygen abundance for the 
metallicity sample. One sigma errors are shown for the oxygen 
abundances and the q parameters. Only the 25 objects with greater 
than one sigma detections of the [OIIIJA4363 line are shown. The 
red line shows the median value of the q parameter in the sample. 



well known it is unfortunately multivalued with both a 
low metallicity and a high metallicity branch. More- 
over, while for the present galaxies we may be reason- 
ably secure that we are on the low metallicity branch 




Fig. 10. — The oxygen abundance versus the R23 parameter for 
the [O III] selected sample (black squares) and the Ha selected ob- 
jects (blue diamonds). Galaxies with [OIII]A4363 lines below the 
two sigma level are shown with open squares at a nominal O abun- 
dance. One sigma errors are shown for the oxygen abundances. 
A linear fit to the combined data sets is sh own by the black line 
while fits to the local data bv |Naeao et al.l |[2006D are shown with 
the green dashed line and by Yin et al.' (200^) as the red dotted 
line. Local data from Izotov and Thuan. (.20071 ) are shown by the 
red triangles. 



(12 + log(0/H) < 8.4), R23 is only weakly dependent 
on metallicity on this branch and has a strong ionization 
dependence (e.g., [McGaugh 1991). 

Nevertheless, recent analyses of substantial samples 
of local galaxies with well determined abundances 
from the direct method have shown a good empiri- 
cal correlation of 12-|-log(0/H ) with R 23 in t he low 
metalhcity r ange (|Nagao et ahl [2006; Y in et al.l [2007). 
Ilzotov and Thuan (2007") favor the Yin et al.l parameter- 
ization 12-hlog(0/H) = 6.486-hl.401x log(R23) based on 
comparisons with their local sample. However, in com- 
paring these to the present sample we must be concerned 
about differences in the galaxy properties and in particu- 
lar whether evolution in the distribution of the ionization 
parameter between the local and distant samples might 
change the relation. 

In Figure [TO] we show the dependence of the present 
measurements of 12+log(0/H) on th e R23 parameter. 
We al so show the local points from Izotov and Thuanl 
(|2007l ) and the empirical fits of Yin et al. {20Q^) (red dot- 
ted curve) and lNagao et al.l (|2606h (green dashed curve). 
It is clear that the present data have lower 12-|-log(0/H) 
determinations at the same R23. This would be expected 
if the ionization parameters are higher in the present 
sample. A factor of five increase in the q parameter could 
easily produce the offsets seen. Empirically we find a re- 
lation of the form 12-hlog(0/H)=6.45 + 0.15 R23, which 
is shown as the black line in Figure [TOl provides a rea- 
sonable description of the data. 

5.3. Metallicity versus H(3 equivalent width 

As is well known the H/3 equivalent width can give 
a rough estimate of the age of the star formation in a 
galaxy. For a Salpeter IIMF and a constant star forma- 
tion rate, EW(H/3) would drop smoothly to a value of 
30 A at about 10^ yr ([Leitherer et al.lll999D while an in- 
stantaneous starburst would drop below this value after 
about 10*" yr. It is therefore of considerable interest to de- 
termine the relation between 12-|-log(0/II) and EW(II/3) 
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REST EW A (IMAGING) 

Fig. 11. — Comparison of rest frame equivalent widths deter- 
mined from the imaging data with those determined from the spec- 
tra. Black squares show [OIII]A5007 equivalent widths and purple 
diamonds Ha equivalent widths. The red line shows the expected 
relation. Only objects where the line wavelength would correspond 
to a transmission above 80% of the peak filter response in the nar- 
row band filter and where the [OIIIJA5007 line is detected in the 
spectra at a signal to noise above 25 are shown. 



in order to find if the metal build up is a function of the 
age of the galaxy. 

We can determine EW(II/3) in two ways. We can mea- 
sure it directly from the spectra or we can determine the 
EW of the line falling in the narrow band filter from the 
imaging observations and then use the line ratios to de- 
termine EW(II/3). However, both methods are somewhat 
problematic. Because the continuua in the spectra are 
extremely weak it can be difficult to precisely measure 
them. The errors here are complex and systemic efforts 
can be important. Thus EW(H/3) can only be confidently 
measured in the highest S/N spectra. The imaging data 
makes a much deeper and more accurate measurement of 
the continuum but translating this to the observed frame 
EW of the measured line requires us to know the filter 
profile extremely accurately. This is a particular prob- 
lem for objects whose wavelengths lie on the edges of the 
filter where the response changes rapidly. We measured 
EW(II/3) with both methods. For the spectral determi- 
nation we restricted ourselves to spectra with S/N> 25 
in the [OIIIJA5007 line. For the images we restricted 
ourselves to galaxies where the transmission in the fil- 
ter was above 80% of the peak transmission. We took 
the continuum from line free broad band colors near the 
narrow band (the / band for the NB816 filter and the z' 
band for the NB912 filters) assuming that the spectrum 
]y was fiat. We determined the EW of the line in the 
narrow band as (lO"'' '*^^"'-^' — where iV is the nar- 
row band magnitude, C is the continuum magnitude and 
</) is the filter response normalized such that the integral 
over the wavelength is unity. Where the line in the fil- 
ter was [OIII]A5007 we converted to EW(II/3) using the 
observed line ratios in the spectra and where the line in 
the filter was Ha we converted to EW(H/3) assuming the 
Case B Balmer ratios. 

We compare the imaging equivalent widths with the 
spectroscopic equivalent widths in Figure 111! There 
is broad overall agreement but a considerable amount 
of scatter reflecting the uncertainties in the two meth- 



ods. However, as is shown in Figure fT2l irrespective of 
which method we use we see a clear relationship between 
12-hlog(0/H) and EW(H/3). The best flt relationship of 
12-hlog(0/H) = 9.72 - 1.07 log(EW(H/3) based on the 
imaging equivalent widths also provides a good descrip- 
tion for the spectroscopically based equivalent widths. 
This relation was less clear in the equival ent widths de- 
termin ed from the poorer spectra used in iKakazu et al.l 
(l200l . With the new data there is now clear evidence 
for the build-up of metals as a function of age in the 
galaxies. 

5.4. Metallicity Luminosity relation 

We next computed the absolute rest frame B mag- 
nitudes using magnitudes from imaging observations in 
bandpasses which are clear of the emission lines and as- 
suming a flat fi, spectral energy distribution to com- 
pute the K correction. We plot these absolute rest 
frame B magnitudes versus the oxygen abundance de- 
rived by the direct method in Figure fT3l The mctallicity- 
luminosity relation can be written as 12+log(0/H) = 
8.01 - 0.15 X {Mb{AB) - 20). 

We can compare the metallicity luminosity rela- 
tion with the l o cal m etallicity luminosity relation from 
iTremonti et"all (|2004') (the green line in Figure I13p and 
with observat ions of magnitude l imited samples at z — 
0.6 — 0.9 from lCowie and Bargerj ()2008[ ) (red points) and 
the corresponding metal-luminosity relation at this red- 
shift (red line). The present sample lies about 0.6 dex 
below the z = 0.6 — 0.9 relation and about 0.8 dex lower 
than the local relation. The Cowie and Barger metal- 
licities are based on the upper branch of the O versus 
R23 relation and may miss some low metallicity objects. 
We have also searched their magnitude limited sample 
for strong emitters (the blue diamonds in Figure [T3| and 
checked these galaxies for [O HI] A4363 but none of these 
show the auroral line. 

6. LUMINOSITY FUNCTIONS AND THE LOCAL 
LYMAN ALPHA EMITTER POPULATION 

We next constructed the Ha luminosity functions 
of the USEL sampl e following the procedures used in 
iKakazu et al.l (|2007f ) but using the present larger spec- 
troscopic sample. We used the narrow band magnitudes 
to determine the line strength of the selection line; Ha 
at low redshift and [O III]A5007 at high rcdshift. For the 
higher redshift objects we computed the H/3 flux from the 
[OIII]A5007 fiux and converted this to an Ha flux using 
the Case B ratio. Because of the high observed frame 
equivalent widths the primary fluxes are insensitive to 
the continuum determination. However, they do depend 
on the fllter response at the emission line wavelength so 
we restricted ourselves to redshifts where the nominal fll- 
ter response is greater than 80% of the peak value follow- 
ing the procedure used in the equivalent width analysis. 
This also has the advantage of providing a uniform se- 
lection and we assume the window function is flat over 
the defined redshift range. Now the volume is simply de- 
fined by the selected redshift range for all objects above 
the minimum luminosity which we take as correspond- 
ing to an observed flux of 1.5 x 10""'^^ erg cm~^ in 
the Ha line. For the high rcdshift sample selected in the 
[O III] A5007 line the true flux limit will generally be lower 
than this since the [OIII]A5007 line is normally stronger 
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Fig. 12. — (a) 12+log(0/H) versus rest frame H/3 equivalent 
widths determined from the imaging data. Black squares show 
[OIII]A5007 selected objects and purple diamonds Ha selected 
objects. The red line shows the best fit linear relation between 
12+log(0/H) and log(EW(H/3)). Only objects with filter trans- 
mission above 80% of the peak and with signal to noise above 15 
for the H/3 line are shown, (b) Same as (a) but with equivalent 
widths determined only from the spectra. All objects with signal 
to noise above 20 for the H/3 line are shown. The red line is the fit 
to the imaging equivalent widths shown in (a). 
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Fig. 14. — The luminosity function of H« at 2 = 0.3 (top panel 
a) and at 2 = 0.7 (bottom panel b). In each case the open symbols 
show the luminosity functions determined from the spectroscopic 
sample alone while the solid symbols show the function corrected 
for the incompleteness in the spectroscopic identification. The er- 
rors are plus and minus 1 sigma and at the highest luminosity we 
show the 1 sigma upper limit. In panel (a) we also show the z = 0.3 
Ha luminosity function from lTresse He Maddoxl l |1998l ) as the blue 
dashed line. The dotted blue line shows the Tresse and Mad- 
dox function multiplied by 0.04 which approximately matches the 
USEL Ho luminosity function. In panel (b) we show the Tresse and 
Maddox function multiplied by 0.18 which approximately matches 
the USEL Ha luminosity function at the higher redshift. 
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Fig. 13. — The oxygen abundance versus the absolute rest frame 
B magnitude for the [O III] selected sample (black squares) The 
red diamonds show the oxygen abundances for a magnitude lim- 
ited s ample in the same redshift range from ICowie and Barged 
1I2OO8I). The green li ne shows the local metal-luminosity relation 
from ITremonti et al.l (12 004'). the red th e metallicity relation at 
2 = 0.8 fronT ICowiean d Bargcij 1 )20081) . and the black line the 
metallicity-luminosity relation of the present sample. 



than Ha and we restrict the sample to objects with Ha 
hne fluxes lying above this threshold. The luminosity 
function is now obtained by dividing the number of ob- 
jects in each luminosity bin by the volume. The incom- 
pleteness corrected luminosity function is obtained from 
the sum of the weights in each luminosity bin divided by 
the volume. Here the weight of an object of a given mag- 
nitude correponds to the total number of objects at that 
magnitude divided by the number of identified ojects at 
that magnitude. Because of the high spectroscopic com- 
pleteness the incompleteness corrections are small except 
at the very lowest luminosities. The 1 sigma errors shown 
are calculated from the Poissonian errors based on the 
number of spectroscopically identified objects in the bin. 
The calculated Ha luminosity function is shown for the 
z = 0.2 — 0.45 range corresponding to the NB816 and 
NB912 Ha selections in Figure fMl a) and the correspond- 
ing Ha luminosity function at z = 0.6 — 0.9 from the 
[OIII]A5007 selected samples in Figure [WbV 

The USEL Ha function at z = 0.3 corresponds to 
abou t 4% of the total Ha lum inosity function at this red- 
shift (jTresse fc Maddol3ll998[ ) which is shown as the blue 
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Fig. 15. — The hya luminosity functions computed from the 
present samples assuming f(Lya)/f(H«)= 5 is shown at 2 = 0.3 
(black squares) and at z = 0.7 (blue triangles). Only objects with 
rest frame equivalent widths above 20A are included. The error 
bars are one sigma computed from the Poisson errors corresponding 
to the number of spectroscopically observed objects in the bin. We 
also show the one sigma upper limit (downward pointing arrow) at 
the highest luminosity. The hya luminosity functi on at the z = 
0.3 re dshift derived from GALEX observations by De harveng et all 
(|2008l ^ is shown with the open red diamonds (actual obervations) 
and solid red diamonds (incompleteness corrected). It is clear that 
there are few local high luminosity LAEs in an area of the present 
size and there is no strong overlap between the present sample and 
the GALEX sample. We also show the z = 3 LAE luminosity 
function of Gronwall et al. (2007) (solid red line) and the same 
function with a downward scaling of 4 (dashed red line) and 20 
(dotted red line) in the number density. 

dashed line in Figure fMl a) . This is similar to the fraction 
of the total star formation rate that Kaka zu et al.l (j2007t ) 
estimated was in the USELs. The USEL Ha luminosity 
function rises rapidly with z = 0.3 paralleling the rapid 
rise in the star formation rates over this redshift interval. 

These Ha luminosity functions also allow us to con- 
struct upper limits on the z — — 1 Lyman alpha emitter 
luminosity functions for comparson wi th the local LAE 
lumin osity function from GALEX (jPeharveng et al] 
|2008( ) and with the measurements of the LAE luminosity 
function at high redshift (z = 2 — 7). As we discuss in 
more detail below a galaxy with a high equivalent width 
in the Lya line will also have a high equivalent width 
in Ha, so the LAEs are a subset of the present sample. 
This allows us to construct a upper bound on the com- 
plete LAE sample to a fixed Lya equivalent width which 
can then be compared to the local and high redshift pop- 
ulation in detail. 

As is well known, the short path length to optical scat- 
tering on neutral hydrogen ensures that Lyman alpha 
photons follow a complex escape from a galaxy. While 
the low metallicity and extinction of the present USEL 
sample are clearly positive indicators of a higher Lya 
escape fraction they are by no means the only factor. In- 
ternal structure, kinematics, and the distribution of star 
formation sites may also play key roles and may in fact 
be the more critical determinants. We take the Case B 
ratio of 8 — 12 for the Lya/Ha flux, depending on the 
elec tron density, to be an upper bound on the line ra- 
tio ([Ferland and Osterbrocklll985( ) (though we note that 
it is possible to have geometries in which the scatter- 
ing can enhan ce the Lya lin e relative to the UV contin- 
uum (Finkelstein et al.ll2007f ).) Given the roughly flat fi, 
continua in these galaxies (or f\ ~ and the rest- 



frame fully complete to a rest equivalent-width of 30A 
compared to the 20A equivalent width normally used to 
select the high-redshift LAE population. The higher red- 
shift z 0.7 sample will contain all LAEs above 20A for 
the case B assumption. However, empirical estimates of 
the Lya escape fraction in the high redshift LAE popula- 
tion made by comparing star formation rates estimated 
from the UV continuum with thos e from the Lya line 
suggest escape fractions of 30-50% (jGawiser et al.l [20071 : 
lGawiserl2008[ ) . If this lower ratio is applicable for the ra- 
tio of the Lya/Ha fluxes, then all LAEs will be included 
in the present samples, though objects close to the Ha 
limiting equivalent width will fall from the LAE sample. 
The LAE sample is therefore a subset of the present sam- 
ple and we can directly obtain an upper bound on the 
LAE luminosity function from the USEL Ha luminosity 
function. 

The sample is small but it provides a first-cut estimate 
of the evolution of the LAE luminosity function. We il- 
lustrate this in Figure 15 where we show the z = 0.3 
and the z = 0.7 LAE luminosity functions that would 
be derived from our Ha sample under the assumption 
Lya/Ha= 5 for objects with rest frame Lya equivalent 
widths above 20A. The z — 0.3 LAE luminosity func- 
tion cam be compared with the direct determinations of 
the LAE luminosity fun ction at this redshift fr om the 
GALEX observations of iDeharveng et all (|2008f ) which 
are shown with the open red diamonds. The GALEX 
data only contain objects with strong UV continuum 
detections and these correspond to the highest lumi- 
nosity emitters. In this sense they parallel the high 
redshift Lyman break galaxies with strong Lyman al- 
pha emission rath er than the line selected LAE samples. 
IDeharveng et al.l estimated that there should be a large 
incompleteness correction to allow for this selection ef- 
fect and the solid red diamonds show their incomplete- 
ness corrected luminosity function. However, the present 
data, even under the case B assumption, suggest that the 
incompleteness c orrections at the low lum inosity end are 
smaller than the IDeharveng et al.l (j2008l ) estimate. (We 
recall that the present data provide an absolute upper 
limit for the case B assumption.) The present z = 0.3 
sample has too small a volume to probe the high lumi- 
nosity end. 

At high redshifts the LAE luminosity function is sur- 
prisingly invariant in the range z = 2.5 — 6 but there are 
signs it begins to drop at both higher and lower redshift. 
The present data show that it must have dropped sub- 
stantially at redshift z — 0.7 and even further at z — 0.3. 
In Figure [TSj we show the most recent determina tion of 
the z = 3 LAE luminosity function from Gronwall et al.l 
(|2007f ). In order to match the present data the LAE lu- 
minosity function would have to drop by about a factor 
of four between z = 3 and z = 0.7 (dashed red line) 
and by a factor of very roughly twenty between z = 3 
and z — 0.3 (dotted red line). Even under the case B 
assumption these drops would be two and ten respec- 
tively. The drop in the LAE luminosity fuction parallels 
the drop in the star formation over this redshift interval 
and there is even a suggestion of donwsizing in the shape 
of the LAE LF at the lowest redshift which is deficient 
in high luminosity objects relative to the low luminosity 
end. 
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7. SUMMARY 

We have described the results of deep spectroscopic ob- 
servations of a narrow band selected sample of extreme 
emission line objects. The results show that such objects 
are common in the z = — 1 rcdshift interval and that 
a very large fraction of the strong emitters are detected 
in the [O III]A4363 line where oxygen abundances can be 
measured using the direct method. The abundances pri- 
marily lie in the 12+log(0 /H) range of 7—8 characteristic 
of XMPGs. We have determined the mctal-himinosity re- 
lation for this class of object finding it lies about 0.6 dex 
below a magnitude selected sample in the same redshift 
interval. We give an emprical relation between R23 and 
12 + log(0/H) which differs from local estimates. We 
also show that low mctallicity objects can be picked out 
by the weakness of [O II]A3727/[OIII]A5007 in the spec- 
tra. The two lowest mctallicity galaxies in the sample 
have 12 + log(0/H) = 6.97±0.17 and 7.25±0.03, making 



them among the lowest mctallicity galaxies known, but 
we expect that as the sample size is increased yet lower 
mctallicity galaxies may be found and that we may hope 
to be able to determine if there is a floor in the galaxy 
mctallicity at these redshifts. 
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